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Abstract 

The aim of this paper is to present a systematic method for computing moments 
of matrix elements taken from circular orthogonal ensembles (COE). The formula 
is given as a sum of Weingarten functions for orthogonal groups but the technique 
for its proof involves Weingarten calculus for unitary groups. As an application, 
explicit expressions for the moments of a single matrix element of a COE matrix 
are given. 
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1 Introduction and main result 

The study of distributions for matrix elements of random matrices is one of main themes 
in random matrix theory. When a random matrix X = (a^y)i<i,i<JV is given, we would 
like to know how to compute the moments 

EjXjjj^Xj^ • • • %i n j„] Or 'K[Xi 1 j 1 Xi 2 j 2 ■ ■ ■ %i n j„X kih x k2h ' ' ' x k m l m }- 

There is a large body of research in this area, but we shall focus on recent papers: 

• Gaussian matrix: There are three typical Gaussian matrix ensembles: Gaussian 
orthogonal, unitary, and symplectic ensembles (GOE, GUE, and GSE); see, e.g., 
|Meh] . The calculus of moments for these is well known, under the name Wick 
calculus; see, e.g., [Z]. 

• Haar-distributed unitary matrix: The compact Lie group U(N) of N x N unitary 
matrices has the Haar probability measure. Collins [C] showed how to compute the 
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moments and called his technique Weingarten calculus, a reference to the trailblazing 
work of Weingarten [W]. The main tool in this theory is the unitary Weingarten 
function Wg^ N \ which is a class function on the symmetric group S n . See also 
[BlIMNj. 

• Haar-distributed orthogonal matrix: The compact Lie group O(N) of N x N real 
orthogonal matrices also has the Haar probability measure. Weingarten calculus 
for O(N) has been developed in [CSl \CM\ IMatlj . In this calculus, the necessary 
function is the orthogonal Weingarten function Wg^ 7V - ) , which is a function on S^n 
with an invariant property under the two-sided action of the hyperoctahedral group. 

• Central Wishart matrix : A Wishart random matrix is a matrix version of a chi- 
square distribution. Graczyk, Letac, and Massam |GLMll IGLM2j established the 
moment computations for the complex and real cases. 

• Inverse of a central Wishart matrix: The distribution of the inverse W~ x of W is 
quite different from the original Wishart matrix W . The inverse distributions of the 
complex and real cases were studied in [GLMlJ and [Mat2j, respectively. For reasons 
not fully understood, unitary and orthogonal Weingarten functions are necessary for 
the complex and real cases, respectively. 

• Noncentral Wishart matrix: The noncentral version of the Wishart matrix is derived 
from noncentral Gaussian vectors. Kuriki and Numata [KN] showed computation 
of the moments using matchings. Note that the problem for the inverted noncentral 
Wishart matrix is still open. 

There are three much-studied circular ensembles of random unitary matrices, intro- 
duced by Dyson: circular orthogonal, unitary, and symplectic ensembles (COE, CUE, 
and CSE). The density function for their eigenvalues Ai, A 2 , . . . , Ajv is proportional to 
rii<i<j<iv l A * - A if with = 1, 2, 4 for COE, CUE and CSE, respectively. See [Mehj for 
details. 

A CUE matrix, that is, a random matrix taken from the CUE, is simply a Haar- 
distributed unitary matrix from U(N), and such a matrix has been studied as stated 
above. Therefore, the next obvious object of study is the COE. The COE is the probability 
space of symmetric unitary matrices defined by the property of being invariant under 
automorphisms V — > t UoVUo, where Uq is a fixed unitary matrix. The COE is a realization 
of the compact symmetric space U(N)/0(N) equipped with the probability measure 
derived from the Haar measure on U(N). 

Our aim in the present paper is to establish the calculus for the moments of matrix 
elements of COE matrices. If U is an iV x iV Haar-distributed unitary matrix, then an 
N x N COE matrix V can be written as V = t UU. In other words, if V = (v^) and 
U = (uij) then = X^fcLi UkiUkj- Therefore, the desired moment for V can be written as 
a sum of moments for U. Since the computation of moments for U is known, the moments 
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for V can in principle be computed. (In this direction, Collins and Stolz |CStj studied 
asymptotic behaviors for the COE and for random matrices associated with classical 
compact symmetric spaces.) Herein, we will express the COE moment in terms of unitary 
Weingarten functions Wg^^ (Proposition [5]) . However, this expression is ill-suited to 
practical computation, so obtaining this expression is not our main goal in the present 
paper. 

Our main result is given as follows. Put [N] = {1, 2, ...,N} and define the right 
action of the symmetric group S 2n on [iV] 2n by 

* CT = (V(i),V(2), • • • ,V(2n)) for * = (*i>*2, • • -,ian) e [N] 2n and a G S 2n . 
Theorem 1. Let N be a positive integer, and let V = (vij)i<ij<N 

be an N x N COE 

matrix. Let i = (ii, . . . , i 2n ) and j = (ji, . . . , j 2n ) be two sequences in [N] 2n , and let 

^f/v(*J 3) = ^[ V il,i2 V i3M ' ' ' V i2n-l,i2n V jl,j2 V j3,ji ' ' ' V j2n-l,j2n] 

and 



M(i,j;z) = J2 Wg^(a;z) (z G C), 



where the sum runs over all permutations a G S 2n satisfying j = i a , and Wg^(a;z) is 
the orthogonal Weingarten function defined in (12.31) below. Then we have 

M N (i,j) = M(i,j;N+l). 
If N + 1 > n, our formula can be written as 

M N (iJ)= Wg° (JV+1) (<x) 

with Wg° {N+1) (a) : = Wg°(<x;iV + 1). However, if N + 1 < n and a G S 2n , then the 
complex function z i-> Wg^(<r; z) may have a pole at z = N + 1, and so Wg^(a; iV + 1) 
does not make sense. Nevertheless, N + 1) is well defined after some cancellations 

of poles. 

Example 1. Let us compute E[|t>n| 6 ] = ElvuVuVnVuVuVn] = M N (i,i) with i = 
(1, 1, 1, 1, 1, 1). Theorem [T] implies M N (i, i) = [^2 aeS Wg^(cr; ^)]«=jv+i- A representation- 
theoretic discussion (the proof of Theorem [2]) gives ^2 a£S(j Wg 3 (a; z) = z ^ z+2 ^ z+ ^ j hence, 

6 

E[|i;i11 ] = (N + 1)(N + 3)(N + 5)' 

If AT + 1 < n = 3, i.e., if AT = 1, we cannot replace E ctGS Wg^(cx; 2:)] z= at + i =2 by 
£ CT£ s 6 W g?(^ 2 ) since Wg^(id 6 ; 2) = z{z+2){ z z ^ z { ;\ z _ 2) has a pole at z = 2. 



4 



Sho MATSUMOTO 



Surprisingly, the orthogonal Weingarten function has appeared three times now in 
random matrix theory: the first appearance was, just as its name implies, in the study of 
Haar-distributed orthogonal matrices [CM] (implicitly |CS] ) and the second was related 
to inverted central real Wishart matrices [Mat 2] . Furthermore, the reader should note 
that the function Wg^ N+1 ' appears, not Wg^^, in our formula. The formula is derived 
by using Weingarten calculus for U(N) and harmonic analysis for symmetric groups, and 
the function Wg° (7V+1) then appears as a result. Since the COE is identified with the 
compact symmetric space U(N)/0(N), it is small wonder that an orthogonal Weingarten 
function appears. However, the theoretical reason why we need 0(N + 1) rather than 
0(N) seems quite mysterious. 

The present paper is organized as follows. In section [2J we review the definition of 
orthogonal Weingarten functions developed in |CMl IMatl[ IMat2j . In section [31 we prove 
Theorem [TJ In section [H using Theorem [Tj, we derive explicit expressions for the moments 
E[|-Ujj| 2n ] of a single matrix element Vij. In section[5j we show that there exist nonnegative 
integers a and b satisfying 

M N (i,j) = aN- n - bN-" 1 - 1 + 0{N- n - 2 ), N -> oo, 

and evaluate a, b. In the final short section, section [6j we show how to use Theorem [I] for 
cases where the degree of integrand n is small. 

Theorem [1] is a very powerful tool for computing various averages on the COE in the 
same way that Wick calculus for Gaussian ensembles and Weingarten calculus for classical 
groups are. We can expect many applications of Theorem [JJ to appear in future research. 

Remark 1. In a recent preprint [Mat 3] . the author showed obtaining some results for 
simple cases of M N (i,j) by the same method as in the present paper. Results in [Mat3j 
seem unsatisfactory compared to the present ones. All these previous results are included 
in the present note with much simpler proofs. 

Remark 2. The COE is a random matrix ensemble associated with the compact symmet- 
ric space U(N)/0(N), which is of Cartan class AI. Other random matrices from classical 
compact symmetric spaces are expressed in terms of Haar-distributed matrices from clas- 
sical Lie groups. Thus, if we apply Weingarten calculus for classical groups to these, we 
obtain the same formulas as those given in Theorem [TJ We leave them as future research. 

Acknowledgment. The author would like to acknowledge an inspiring conversations with 
Benoit Collins. This work was supported by a Grant-in- Aid for Young Scientists (B) No. 
22740060. 

2 Orthogonal Weingarten functions 

We first review the theory of Weingarten functions for orthogonal groups; see |CM[ IMatl|, 
IMat2] for details. Claims in the first two subsections can be found in [Mac[ VII. 2]. 
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2.1 Hyperoctahedral groups and coset-types 

Let H n be the subgroup in SVn generated by transpositions (2k — 1 2k), 1 < k < n, and 
by double transpositions (2i — 1 2j — l)(2i 2j), 1 < i < j < n. The group H n is called 
the hyperoctahedral group. Note that \H n \ = 2 n n\, and the pair (S2 n ,H n ) is a Gelfand 
pair in the sense of [Macl VII]. 

A partition A = (Ai, A 2 , . . . ) is a weakly decreasing sequence of nonnegative integers 
such that |A| := Yli>i * s finite. If |A| = n, we call A a partition of n and write A h n. 
The length (( A) of A is defined as the number of nonzero Aj. 

Given o G Sim we attach an undirected graph r(<r) with vertices 1, 2, . . . , 2n and edge 
set consisting of |{2fc — 1,2k} \ k G [n]} and |{cr(2fc — l),cr(2A;)} | k G [«-]}• Each vertex 
of the graph lies on exactly two edges, and the number of vertices in each connected 
component is even. If the numbers of vertices are 2Ai > 2A 2 > • • • > 2A/, then the 
sequence A = (Ai, A 2 , . . . , \i) is a partition of n. We will refer to A as the coset-type of cr. 

In general, given cr, r G S^ni their coset-types coincide if and only if H n crH n = H n rH n . 
Hence, we have the double coset decomposition of H n in 5 , 2n 

5* 2n = | | H^, where = {cr G | the coset-type of a is /i}. 

fi\-n 

Note that = and = (2 n n!) 2 /(2^z M ), where 

r>l 

with multiplicities m r (fi) — \{i > 1 \ fa — r}\. 

Denote by (-'(cr) the number of connected components of r(cr), or equivalently, ('(cr) = 
e(ji) if c G Hp. 

2.2 Zonal spherical functions 

For two complex- valued functions fi, / 2 on Sin, their convolution /i * / 2 is defined by 

(fl * / 2 )(*) = X) /iK'^W ^ G S *0- 

Let e# n be the function on defined by 

e H (a) = { (2 " n!)_1 if(TGif "' 
[O otherwise. 

For each A h n, we define the zonal spherical function u x by w A = y 2X * en n = &H n * X 2X , 
or by 

w V) = air E * 2A ^ ( ff G 5 ^)> 
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where \ 2X is the irreducible character of S 2n associated with 2A = (2Aj, 2X 2 , . . . ). The 
functions u x are if n -biinvariant, and so they take a constant value on each H^. 
We will use the following two lemmas later. 

Lemma 1. If X h n and /i h 2n, £/ien cj a *x m = y2jrS 2 \ >fl u) x ■ Here, f 2X is the dimension of 
the irreducible representation of character \ 2X ■ Hence, u x * u p = y^-5 Xp uj x for X,p\~n. 

Proof. The well-known fact x M * X p = ^jjr^p,,pX p gives cu x * \ p = e H n * X 2X * X p = 

(2n)! r 2A (2n)! r X n 

Lemma 2. For any complex number z and a G S2 n , 

(2.1) ^ = SE^yw 

^ Ahn 

(2.2) c x (z)= n^^-'-^nn^^-^ 1 )' 

(i,j)e\ i=l 3=1 

Proof. This is a special case of an identity involving power-sum symmetric polynomials 
and zonal polynomials. See [Mat2l (4.9)]. □ 

2.3 Orthogonal Weingarten functions 

For each positive integer n, we define the orthogonal Weingarten function as 

2 n n' f 2X 

(2.3) Wg>; z) = -—Y. ^77f V) (* G S 2n , z G C), 

^ Ahn A ^ ' 

which was first defined in |CM] . If we fix <r, the complex function z Wg^(cr; z) is a 
rational function in z. Since C^(z) has a factor (z— — 2) ■ ■ ■ (z — £(A) + 1), the function 
Wg^(a; z) may have poles at the points z — 1, 2, . . . , n — 1 on the half line R >0 . 
If JV > n, then C' X (N) ^ for all Ahn, and we can define 



2A 



?2nJ 



Wg?"V) = [Wg?(a; z)) z=N = ^E cpO wA(<r) (a G 5: 

We will refer to the above as the Weingarten function for the orthogonal group O(N). 

The function S 2n 3 o ^ Wg^(cr;z) is constant on each double coset H^. We will 
denote by Wg n (fi; z) the value of Wg^(o~; z) at a G -f^. 
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Example 2. The following examples can be found in |CM[ ICS] . 

Wg?([l];s) = i 

Z 



Wg 2 °([2]; 2) = -— -- -, Wg?([l 2 ]; z 



+ 2)0-1)' J ' y z(z + 2)0-1) 



Wg3( -'- z ( z + 2)0 + 4)0 -1)0 -2)' 
Wg?([2, !];*)= ' 



^0 + 4)0- 1)0-2) ! 
Wg?([l 3 ];,) : 2 + 3:-2 



^0 + 2)0 + 4)0- 1)0-2)' 



Remark 3. The orthogonal Weingarten function was introduced in [CS] implicitly, in a 
study of moments for Haar-distributed orthogonal matrices. The explicit definition was 
first given in |CMj . If O = 0ij)i<M<iv is a Haar-distributed orthogonal matrix taken from 
O(N) and N > n, then for each a G S 2n , 

Wg^ (A °0) = "l ■ i.' <> - .i.. <> -:i ; • • • °n,n n -i°n,hn] 

with (ji,j 2 ,---,j2n) = (l,l,2,2...,n,n) CT . 



3 Proof of Theorem [T] 

3.1 Haar-distributed unitary matrix 

The orthogonal Weingarten function appears in the statement of Theorem [H but we must 
employ a unitary Weingarten function in our proof. 
We define the unitary Weingarten function as 

(3.1) Wg^O; z) = - Y, tXtxV) {° eS n , zeC). 

Here, x X is the irreducible character of S n associated with A, f x is its degree, and 

As in the case of the orthogonal Weingarten function, the function Wg^(cr; z) may have 
poles at the points z = 1, 2, . . . , n — 1 on R>o- 
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If N > n, we let 

Wg?^(a) = [WgJ(a;«)W 



1 E 



m ~ C X (N) 



and refer to it as the Weingarten function for the unitary group U(N). 



Given a unitary matrix U = (u 
= (ji, • • • ,jn), *' 



(*!)•••) *m)' 



) G £7(iV) and four sequences i = 
(ji) - " ifm) °f indices in [iV], we let 



^i,i 2 , • • ■ ,i n ), 



U(i,j\i,f) = ", j J',,,, ■ ■■Ui njn u l[j/ u i ^ ■ •■u i!mj , n . 

Lemma 3 (Weingarten formula for the unitary group; see, e.g., [UJ ICSj ). Let U be 
an N x N Haar- distributed unitary matrix and let i = {i\,i2, ■ ■ ■ ,in)> 3 = (ji? • • • ->3n)> 
i' — (ij, . . . , i' m ) , j' = (ji, . . . , j' m ) be four sequences of indices in [N] . If n = m, 



E[U(i,j\i',j')} 



E E W gn(^V;,) 



i"=i'j T =j' 



z=N 



and it vanishes otherwise. 

Example 3. If U = (iHj) is an N x iV Haar-distributed unitary matrix, then 

[2WgJ(id 2 ;z) + 2Wg?((12);*)] 



E[\u u \ 
2- — 



E E w « 

(T6S2 T&S2 
1 -1 

T7 7T + 2- — 



=7V 



z=N 







r 2 " 


2 


1)- 


z=7V 


Vz{z + l). 


2 =7V ~~ N(N + 1) 



for all N > 1. 

3.2 Expression in terms of unitary Weingarten functions 

Let V = {vij)i<i ) j<N be an N x N COE matrix. Let M(n) = {Si )ir (j))i<i,j<N be the 
permutation matrix associated to 7r 6 SV- Since t M(7r)VM(n) has the same distribution 
as V, we have 



V12 (i^j) 



where X = Y means that X has the same distribution as Y. 

Using a Haar-distributed unitary matrix U = (t%-)i<ij<iv, we can write V = t UU. Let 
i = (ii,i 2 , ■ ■ .,h n ) e [N] 2n and j = ■ ■ ■ , hw) € [N] 2m . Since Vij = J2k u ki u kj, we 

have 



(3.2) E[t) ilii2 t) i3i 



•Vi 



»2tt-li*2n V hd2 V 33,34, "Jlm-ldlm 



E 

fee [JV]« «e[iv] m 



E E[C/(fe,i|i,j)], 
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where 



k — {k\, k2, • • • , k n ), 

k =(kij kx, k 2 , k2, ■ ■ ■ , k n , k n ), 

Using f !3.2p and Lemma El gives the following. 



I —{h, h, ■ ■ ■ , lm), 

I =(h, h, h, h, ■ ■ ■ , lm, lr, 



Lemma 4. E[v iui2 v i3ii4 ■ ■ ■ v^^^v^ j 2 v j3 j 4 ■ ■ ■ v j2m _ uj2m \ vanishes unless n = m and 
there exists a permutation a G S 2n such that j = i a . 

Hereafter, we will assume that i = (ii, . . . , i 2n ) and j = (ji, . . . , j 2n ) and let 
M N (i,j) = EK,i 2 t> 



«3,«4 ' ' ' V i2n-l,i2n V jl,j2 V j3,j4 ' ' ' V j2n-l,j2n\- 



Our first step in the proof of Theorem [T] is to prove the following. 
Proposition 5. For a complex number z, let 

(3.3) M(i,j;z)= £ E ^ (T) Wg 2 ^(r-V;z), 

j = i° 

where £'(t) is the length of the coset-type ofr. Then M N {i,j) = [M(i,j;z)] z=N . 
Proof. Equation (13. 2p . Lemma [31 and a change of the order of sums give 



k£[N] n le[N]" 



u 'r- l a:z) 



J2 J2 w &2, ; >" "■ 

l=(k) T 3=i a 



-N 



E E W < 



E E 1 

ke[N] n l£[N] n J Z =JV 
l=(k) T 



and so the claim follows from Lemma M below. 



□ 



Lemma 6. Let t G S 2n and let A(r,N) = {(k,l) G [N] n x [N) n | I = (k) T }. Then 
\A(t,N)\ = N 1 '^. 

Proof. Given (k,l) G A(t,N), write k = (k\, k 2 , ■ ■ ■ , k n ), I = (1%, l 2 , ■ ■ ■ , l n ), and k = 
(ki, k 2 , k 3 , /c 4 , . . . , k 2n _i, k 2n ). Then it follows immediately that 



(3.4) lj = fc T (2j-i) = k T (2j) and kj = k 2j -i = k 2j 
In particular, I is determined by k and r. 



(j = l,2,...,n). 
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Recall the definition of the graph T(r) associated with r. Given a sequence k e [N] n , 
we assign numbers hi, k 2 , . . . , k 2n to vertices 1, 2, . . . , 2n of T(r). Then the condition (j3.4p 
implies that two vertices r and s belong to the same connected component of T(r) if and 
only if k r = k s . 

For example, let n — 4 and t = (41652873) e ^- In this case, ( 13. 4 p gives the 
chain of two equalities k± = k\ = k 2 = k s = k 7 = k 3 = k 4 , k§ = k 5 = kg, so we have 
\A(t,N)\ = \{ik 1 ,k 2 ,k 3 ,k 4 ) E [iV] 4 I ki = h = k 4 }\ = N 2 . On the other hand, the 
graph T(r) consists of two connected components, one of which has vertices numbered by 
4, 1, 2, 8, 7, 3 and the other has vertices numbered by 6,5. 

Hence, in general, \A(t, N)\ coincides with N t<yT \ where i'{r) is by definition the 
number of connected components of r(r). □ 

3.3 Expression in terms of orthogonal Weingarten functions 

Substituting (EEJ and §n§ into Q gives 

1 v- / M 



,(2n). 

0-e5 2 n T&S 2 n Ahn 



E - t -^- 1 < 



( .,(.:. 1 

fih2n 



(2n)! C„(*) 



2"n! 2X C' x (z) x 

2^ 2^2^ {(2n)\y J 1 c^zy x ){ h 



Tg52n Ahn uh2n 

so Lemma [1] gives 

cr€b2n Ahn 
i=i CT 

Since C 2A (z) = C' x {z)C' x {z + 1), it follows from Q that 

on„| x2A 

= <50f E Ec^TT)-»= E wg?(.;, + i). 

V ' a£S 2n Ahn AV ' <rGS 2n 

Hence, the last equation together with Proposition [5] implies Theorem [TJ 

4 Moments of a single matrix element 

We next obtain explicit values for the moment of a single matrix element E[|v 

Theorem 2. Let N > 1 and let Va be a diagonal entry of an N x N COE matrix. Then, 
for n > 1 , 

2 n nl 

E[K| 2nl - 



|2nl 



(N+l)(N + 3)---(N + 2n- 1) 
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Proof. Theorem [T] implies 



E[\vuH = M N ((i,...,i)M,...,i)) 



2n 



In 



E w s>; 



0-£S2n 



=AT+1 



Since x (2n) = 1, Lemma □ implies that Eaes 2n W V) = (u x * x {2n) )(id 2n ) = (2n)\5 x ,( n ), 
and so we have 



E w ^°^> 



_ 2"n! v / 2A v 



2 n n\f {2n) 



Ahn 

2 n n! 



q n) (z) «(z + 2)...(^ + 2(n-l))- 
which implies the desired result. 



□ 



Theorem 3. Let N > 2 and let &e an off- diagonal entry of an N x N COE matrix. 
Then, for n > 1, 



ft! 



(iv+2n-i)n:: 2 (^+*o 



The proof of Theorem [3] is more difficult than that of Theorem [2j We first show the 
following lemma. 

Lemma 7. Under the same assumptions as in Theorem^, 



(4.1) 



ENu 



1 2m 



71! 



2 _ p2\ ^a( 2 ) 
(2")! ^ ^ + 1)' 



Ahn A ^ 
£(A)<2 



Proof. Let i = . . . ,i,j) G [iV] 2ri , where i ^ j- If ex G satisfies i CT = i, then it 

permutates odd numbers and even numbers, and so it is uniquely expressed as 

a(2k - 1) = 27Ti(Jfe) - 1, a (2k) = 2n 2 (k), k = l,2,...,n 

for some m, tt 2 G S n . In this case, it is not difficult to see that the coset-type of a coincides 
with the cycle-type of ttitt^ 1 - It follows from Theorem [1] that E[|i> ij\ 2n ] = i; z)] z= n + i 

with 

M(i,i;z)=J2 W ^;z)= E Wg°(Cycle(7r 1 7r 2 - 1 );z) 
=n! E W Sn (Cycle (vr);z), 

7TGS„ 
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where Cycle(7r) is the cycle-type of it and Wg^(/x; z) is the value of Wg^(a; z) at a 6 i?^. 
The number of permutations in of cycle-type is n!/z„; hence, 



A4(w)=n!^-Wg°(/z^). 

I 

Since = (2"n!)7(2^z M ), we have 

M(i,i^) = 2- 4 "^2^)|^|Wg«( yU ;z) = 2- 4 " £ 2^Wg>;* 

by the decomposition = LL-n-^/x- Using (12. ip . (12. 3p . and Lemma [T] gives 

V ^ ' ' ' Ahn /ihn A* ^ ' 

_("!) J V^( 2 ) 



where we used a/(id 2n ) = 1 to obtain the last equality. Since C' x {2) = if £(A) > 2, the 
range of the last sum can be restricted to {A h n | £(X) < 2}. Since C' X (N + 1) is nonzero 
for every A in the restricted set, we have 



(2n)! ^ 7 Ctftf + l)" 



E[| % H = [M(M;*)W +1 = ^ £ / 2A ^ 

Ahn 
£(A)<2 

□ 



We next evaluate the sum on the right-hand side of (14.11) . 

If fj, is a partition of n, the well-known hook-length formula states that = n\/h(fi), 
where h(n) is the hook-length product 

= II (^i + Vj - i - J ' + !) 

and fjf = /4, • • • ) is the conjugate partition of 

A partition A of n with £(A) < 2 can be expressed as A = (n — r, r) with some 
< r < [n/2], where [n/2] is the integer part of n/2. We now fix n, N and let 

, _ / {2n - 2r ' 2 ^ ( U r ,)(2) _ q„_,,)(2) 



(2n)! ■ q Mr) (JV + 1) fc((2n - 2r, 2r))C ( ' n _ rjr) (iV + 1) 
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for each < r < [n/2]. A direct computation gives 

2n-4r + l ^ 1 ^ 1 

r ~ (2r)!! ■ (2n - 2r + 1)!! JLJL at + 2j - 1 ' JLJL N + 2j' 

Here, (2ife - 1)!! = IlJ=i(2j - 1) and (2ife)!! = nj=i(2j) with 0!! = 1. 
Lemma 8. For each < k < [n/2], 

1 

A + A x + ■ • ■ + A k 



(2k)\\ ■ (2(n - k) — 1)!! 



j n— fc— 1 ., fc— 1 



n — tt— 

J- J- iV + 2?-l + 



X (N + 2n — 1) JLX TV + 2j - 1 JJL N + 2j' 

Proof. We show the lemma by induction on k. Observe that the case at k — holds true: 

1 " 1 
A ° = (2n-l)!! j2 N + 2j - 1 

We next consider the identity at k + 1. From the induction assumption, we have 
A + A + • • • + A k + 

= (2A;)H ■ (2n-2fc- 1)!! ' (N + 2n — 1) ' II 7V + 2j - 1 ' Q N + 2j 



+ 



2n-4k-3 " '" 1 



n — n— 

J- J- iV + 2?-l + 



(2Jfe + 2)!! • (2n — 2/c — 1)!! ^ iV + 2 j — 1 JA iV + 2j ' 
By a reduction, the last equation is 





1 




(2k + 2)(N + 2k) + (2n - Ak - 


-3)(iV + 2n-l) 


(2fc + 2)! 


\(2n-2k- 


1)! 


!( iV + 2n-l)n;r i fc " 1 (iV + 2j- 


i)-n- =0 (^+2j) 




1 




(2n-2fc- l)(iV + 2n 


-2k -3) 


(2fc + 2)! 


\(2n-2k- 
1 


1)! 


!( iV+ 2n-i)n;rf" 1 (iv+2j- 

1 


i)-n-= (A r +2j) 



(2A; + 2)!!(2n -2k- 3)!! (AT + 2n - 1) njlf ~ 2 (N + 2j - 1) • nJ=o(# + 2 ^') ' 
which concludes the desired identity. □ 
Finally, we will prove Theorem |3j 
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Proof of Theorem^ From (14. Xp and the definition of A r , we have 

[n/2] 

E[\v tl \ 2n } = (n\) 2 J2A. 

r=0 

We first suppose that n is even: n = 2m. Using Lemma El we have 

m ^ ^ m— 1 ^ m— 1 ^ 

^ ^ r ~ (2m)!! • (2m - 1)!! (iV + 2n - 1) ' I| jy + 2 j - 1 ' IJ jy + 2 j 

=_ t rr 1 

(2m)! (N + 2n — 1) ^1 # + jfe' 

which implies Theorem |3] The same result holds for the odd n case by a similar calculation. 
Thus, the proof of Theorem [3] is complete. □ 

Remark 4. Let on and Ojj (i ^ j) be two diagonal entries of a Haar-distributed orthogonal 
matrix from O(N). Theorem 5.1 in [CM] gives that 

E[(<* + 0ii n= £ / 2A ^§r 

^(A)<2 

for all n > 1 and N > 2. From the discussion in this section, we have obtained a new 
identity 

E[(o« + %) 2n ] = ^ ^ • 

11 J n! (iV + 2n-2)n^ (^ + ^-l) 

5 Asymptotic behavior 

We now give an asymptotic expansion for Mjy(i, j) for the case iV — > oo. In this section, 
do not confuse Landou's big O notation 0(N~ k ) with the symbol for the orthogonal group 
O(N). 

As we saw in subsection 12.3^ 

Wg n °^(a) = Wg>; ^) = E "V) (- e * 2n ) 

if iV is sufficiently large. We already know the asymptotic behavior of Wg^^ as N — > oo 
as follows. 

Lemma 9 QMatlJ). Fix n and a G S^. iV — > oo, 
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1. If a G H n , then Wg° (A0 (a) = N ~ n + 0{N~ n - 2 ). 

2. If a G i? (2j i»-2), then Wg% {N) (a) = -N~ n - 1 + 0(N~ n - 2 ). 

3. Otherwise, Wg° (A °(<x) = 0(N~ n - 2 ). 

Theorem 4. Fix a positive integer n and fix two sequences i = . . . , 2 2 n) «^ J 
(ji, ■ ■ -,32n)- As N -> 00, 



M*(i, j) = s(i,j)N- n - (s'(i,j) + sfrjWN-"- 1 + 0(N- n - 2 ), 

where 

s(i,j) = \{aeH n \j = i a }\ and s '(iJ) = \{aeH {2A n-2 ) \j = i a }\. 
Proof. Since (N + l)~ n = N~ n - nN^- 1 + 0(N- n ~ 2 ), the last lemma gives 

{N~ n - nN- n ~ l + 0(N~ n ~ 2 ) if a G tf n , 
—jy-"- 1 _)_ 0(N~ n ~ 2 ) if ff 6V 2 )- 

0(jV-"- 2 ) otherwise. 

It follows from Theorem [1] and the double coset decomposition S^n = LLi-n that 
M N (i,j) is 

=s(i,j)(N~ n -nN- 11 - 1 ) - s'ii.j)^ 71 - 1 + 0{N~ n - 2 ) 
as N — > oo. □ 

If s(i,j) > 0, then there exists a ( G such that j = i^. Fixing £, the correspondence 
t : cr I—)- cr£ defines a bijection i : {a G if n | i = i "} — > {r G | j = i r }. Therefore, in 
this case, s(i,j) coincides with the order of the stabilizer subgroup of i in H n . 

Example 4. 1. Let i = (i, i, . . . , i) G [iV] 2n . Then s(t, i) = \H n \ = 2 n n\ and s'(t, i) = 
\H(2, i"- 2 ) I — n ( n — l)2 n n!, and so we have 

E[|^| 2 ] = 2 n n\{N- n - n 2 ^ 11 - 1 ) + 0(N~ n - 2 ). 

2. Let i ^ j and let i = . . . Then s(i, i) = n\ and s'(i,i) = u)^- 

Therefore, 



E[|^| 2 ™] =n! (V 
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3. Let i = (11,12, ■ ■ ■ , «2n) be a sequence with distinct numbers. Then s(i,i) = 1 and 
s'(i, i) = 0; hence, 

E[K i2 v i3i4 ■ ■ ^ 2n _ 1)i2 J 2 ] = N~ n - nN~ n - 1 + 0(N- n - 2 ). 

4. Let i = and j = {i,i,j,j) with % ^ j. Then s(i,j) = and s'(i,j) = 4, 
and so we have E[vf j v~v~] = -4N~ 3 + 0(N~ A ). 

A standard complex Gaussian random variable Z is defined by Z = (Xi + yJ— 1X 2 ) / \/2, 
where Xi and X 2 are independent and standard (real) Gaussian random variables. Then Z 
satisfies ~E[Z n Z ] = 5 nm n\ for any positive integers n and m. Lemma S] implies Eff^i^" 1 ] = 
E[f^t^J m ] = with i 7^ j if n ^ m, and the first two claims in Example S] (or Theorem El 
[3]) imply 

lim (—\ E\\v ii \ 2n ]= Mm N n E\\v ij \ 2n ]=n\. 

Hence, by the Carleman criterion for moment problems (see the proof of Theorem 3.1 in 
[PR] ) . we obtain an algebraic proof of the following result, which was proved by Jiang [J]. 

Corollary 5. Let = (vj-lp), N > 1, be a sequence of COE matrices. Fix positive 

integers i ^ j. As N — > oo, the random variables a/ ~N~/ 2v\^ and yfNv\^ converge weakly 
to a standard complex Gaussian random variable. 

6 Example 

Finally, we give some examples for small n by applying Theorem [TJ Let V = (vy) be an 
N x N COE matrix as usual. We have already obtained explicit values of moments of the 
form E[|fjj| 2ri ] when describing Theorems [2] and [3j and so let us consider other examples. 
Let i = (ii, i 2 , . . . , «2n) and j = (ji, j 2 , . . . , j2n) be sequences in [N] 2n . Let 

s li (i,j) = \{aeH fl \ j = i a }\ 

for each \i h n. Note that s^(i,j) = s(i,j) and S( 2j in-2)(i, j) = s'{i,j) in the notation 
of Theorem HI Then Theorem [T] implies 

M N {i, j) = ^(»,i)Wg^; N + 1), 

which is an identity in N, and iV is regarded as an indeterminate rather than a fixed 
positive integer here. Some explicit values for Wg^(/i; z) can be seen in Example [2] or in 
[CSllCM] . 
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Example 5. (degree 1.) If h,i 2 ,ji,j 2 are positive integers in [N], 



N- 



- if % X = % 2 = jl = j 2 , 



Hv ili2 v jlh } = M N ((ii,i 2 ), {juh)) = { ]vTT if *i ^ z " 2 > ^ ^ and Oi'^} = {ji, j 2 }, 

otherwise. 

Example 6. (degree 2.) Let i ^ j. Put 

T _ rn 2 3 4\ /l 2 3 4\ / 1 2 3 4\ C1 2 3 4U 
n l — Ul 2 3 4M2134Ml 2 4 3 /A2143/J> 

T _ |n 2 3 4\ C1 2 3 4\ / 1 2 3 4\ ( 1 2 3 4 \\ 
L 2 - Ul 3 2 4j)l3 124Ml 3 4 2 M3 142W' 

Using the list of Wg 2 z) given as Example [2j we can compute 

n\v u v n \ 2 } =M N ((i,i,j,j),(i,i,j,j)) = Wg^(a;iV+l) 



o-eLi 

=4-Wg^([l 2 ];iV + l) 



o, ri2 , R7 , ^ 4(iV + 2) 



(AT + 1)N(N + 3) 



%£■ Wj] =M*((1, 2, 1, 2), (1, 1, 2, 2)) = N + !) 

0-GL2 

=4.Wg°([2l;i V + l) = (jv+1) - 4 (jv + 3) . 

Example 7. (degree 3.) Using Theorem [1] and the list of Wg^(/i; z) given as Example EJ 
we have 

o ( ^. AT , ^ 8(iV 2 + 5iV + 2) 



E[|vii«22V33r] = 8 ■ Wg^([l d ]; N + 1) 



(N - 1)N(N + 1)(N + 3)(JV + 5) 



E^^eWiwiF] =Wg° (({ 1 1 1 1 1) ;N + 1) = Wg°([3];iV + 1) 

2 

~ (iV - 1)N(N + l)(iV + 3)(7V + 5) ' 



E[\v 2 u v 22 \ 2 } =16 • Wg°([l 3 ]; N + 1) + 32 ■ Wg°([2, 1]; N + 1) 

16(iV + 4) 



iV(A^+l)(iV + 3)(iV + 5)' 
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